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Abstract. A recent paper [Jokela et ai, arXiv:0806.1491 (2008)] contains a surmise 
about an expectation value in a Coulomb gas which interacts with an additional 
charge £ that sits at a fixed position. Here I demonstrate the validity of the surmised 
expression and extend it to a certain class of higher cumulants. The calculation is based 
on the analogy to statistical averages in the circular unitary ensemble of random-matrix 
theory and exploits properties of orthogonal polynomials on the unit circle. 
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1. Purpose and result 

In a recent paper Jokela, Jarvinen and Keski-Vakkuri studied n-point functions in 
timelike boundary Liouville theory via the analogy to a Coulomb gas on a unit circle 
PQ. In this analogy, N unit charges at position ti interact with additional charges of 
integer value £ a , situated at position r a . To illustrate this technique the authors of [1] 
considered the canonical expectation value 

(•)-|/fi§n|^-e^i 2 nh iT -e^| 2? (-) w 

i=l i<j i 

(where Z is a normalization factor so that (1) = 1) and surmised that 

(Reax) = (T cos(r - *)\ = - J^L. (2) 

In this communication I demonstrate the validity of (j2j), and also compute 
expectation values of the more general quantities 

a n = ex P (iSte* _ r ) ) • ( 3 ) 

il<i2<---<in \ k=l / 

As a result, I find 

<°-) = (- 1 )" ' (N+ lAiyJ> V„ = 0,1,2,...,JV;{>0, (4) 

where (xY y ' = T(x + y)/T(x) is the generalized rising factorial (Pochhammer symbol). 
In particular, the validity of ([2} follows from (jlj) by setting n = 1. 



Statistical correlations in a Coulomb gas with a test charge 



2 



Expression (J3J) will be obtained by relating the generating polynomial 

N 

tp N ,t(\) = j:(a n )(-\) N - n (5) 

n=0 

to a weighted average of the secular polynomial in the circular unitary ensemble 
(CUE). This in turn establishes a relation to the Szego polynomial of a Toeplitz matrix 
composed of binomial coefficients. This calculation sidesteps Jack polynomials and 
generalized Selberg integrals, which can be used to tackle general expectation values in 
multicomponent Coulomb gases [2]. 



2. Reformulation in terms of random matrices 



The CUE is composed oi N x N dimensional unitary matrices U distributed according 
to the Haar measure. Identify tj with the eigenphases of such a matrix. The joint 
probability distribution is then given by [3] 

nWf=i) = ^n|^-e^| 2 , (6) 

i<j 

where z is again a normalization constant. This expression can also be written as the 
product of two Vandermonde determinants det V + det V~ with matrices = e 10 ^" 1-1 ^. 
Furthermore, we can write 

J] |e ir - e" l | 2? = [det(l - Ue~' 1T ) det(l - C/V^]^ (7) 

i 

Finally, the expressions a n in ([3]) arise as the expansion coefficients of the secular 
polynomial 

N 

det(f/e- ir - A) = o n (-A) JV ~ n . (8) 

n=0 

Note that in all these expressions r can be shifted to any fixed value by a uniform shift 
of all tj's, which leaves the unitary ensemble invariant. Therefore the expectation values 
are independent of r. Collecting all results, we have the identity 



Jdet(l-{/)det(l-t/t)| Sd e t ([/-A) /c| 
WA) " <[det(l-f/)det(l-f/t)] { ) CUE ' (9) 

This can be interpreted as a weighted average of the secular polynomial in the CUE. 



3. Random-matrix average 

Statistical properties of the secular polynomial without the weight factor (£ = 0) have 
been considered in pTJ. Clearly, ipN,o = (— A) , so that in this case the attention quickly 
moves on to higher moments of the a n . The main technical observation in [3] which 
allows to address the case of finite £ concerns averages of expressions g({ti\f =l ) that are 
completely symmetric in all eigenphases. In this situation the average can be found via 

r r\t- 

(<?(0*}£Li)>cue = J II ^9({ti}li) det W, (10) 
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where Wi m = e ltm "~ m h Equation (1101) is simpler than the general expression involving 
the product of two Vandermonde matrices, since each eigenphase only appears in a 
single column of W. 

In the present problem, the numerator in (Q is represented by the completely 
symmetric function 

9i({ti}ti) = f[i(d U - A)(l - e*)*(l - e-^], (11) 
i=i 

while for the denominator we need to consider the similar expression 

TV 

92mii) = m i '^( i -^n ( i2 ) 

Using the multilinearity of the determinant we can now pull each factor into the 
ith column and perform the integrals. This delivers the representation 

detCB-AA) 

<P*M = detA , (13) 

where the matrices A lm = (-l) z_m ( e+ jl m ) , B lm = (-l) z ~ m+1 ( e+ ^ +1 ) have entries 
given by binomial coefficients. We now exploit the regular structure of these matrices 
in two steps. 

1) Matrix B contains the same entries as matrix A, but shifted to the left by one 
column index. In order to exploit this, let us expand the determinant in the numerator 
into a sum of determinants of matrices labeled by X = (x m )^ =1 , where we select each 
column either from A (x m = A) or from B (x m = B). [Note that we set these symbols 
in roman letters.] The related structure of A and B then entails that det X vanishes if 
X contains a subsequence (x m ,x m +i) = (A, B). Consequently we only need to consider 
determinants of matrices X n = (B)^ =1 ©(A)^ =n+1 , associated to sequences that contain 
n leading B's and N — n trailing A's. As A is multiplied by —A, det X n contributes to 
order (-A)*"". (Note that X = A and X N = B.) 

2) Next, consider the matrix An + \, where the subscript denotes the dimension, and 
strike out the first row and the n + 1st column (n — 0, 1, 2, • • • , N). This takes exactly 
the form of the matrix X n of dimension N. Therefore, the expressions (— l) n det X n 
are the cofactors of the first row of A^+i- These, in turn, are proportional to the first 
column of A^- +1 , where the proportionality factor is given by det A/v+i- Consequently, 
taking care of all alternating signs, 

det An n=0 

Via steps 1) and 2) we have eliminated any reference to the matrix B. 
4. Orthogonal polynomials 

Matrix A is a Toeplitz matrix, Ai m = ci_ m . In order to find the explicit expression (j3J) we 
now make contact to the theory of orthogonal polynomials on the unit circle [5]. Among 
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its many applications, this theory provides a general expression for the inverse of any 
Toeplitz matrix in terms of Szego polynomials For the case of real symmetric 

coefficients, the inverse is generated via 

X^ipNjX^N^- 1 ) - X N fiip N {X- 1 )ilj N {fi) _ det A N+1 " x N _ n m 

T - i I A K A N+1 )m+l,n+lA V ■ [i-O) 



In 



Comparison of this equation with m = to ()T4j) immediately leads to the identification 

of (— l) N (p N ^(X) with the Szego polynomial Viv(A) of degree N. These polynomials 

satisfy recursion relations which for real symmetric coefficients take the form 

1 r c\X 00 
-jj— f ^n-i(X) £ CnX\ (16a) 
Qn-i J 2tti n fr' 0O 

Vjv(A) = A^at-i(A) + lN X N - 1 i) N _ 1 {X~ l ), (166) 

b~N = S N -i(l - 7^). (16c) 

The initial conditions are So = cq, ipo(X) = 1. The numbers 7^ are known as the Schur 
or Verblunsky coefficients. 

It can now be seen in an explicit if tedious calculation that the polynomials 

*Pn(X) = (-1) ^(A) = X, (iV + l)(g)(l)«-D A (17a) 

= X N 2 F 1 (-N,Z;-N-Z;X- 1 ) (176) 

[with coefficients and expansion given in (TJ|, (j3J)] indeed fulfill the Szego recursion 
generated by the binomial coefficients c n = (— l) n Q^ n ). The recursion coefficients take 
the simple form 

This completes the proof of (j4l, and also entails the validity of (J2]). 
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